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ABSTRACT 

Although calculations of the b — > 57 rate in supersymmetric grand unified 
models have always either ignored the gluino mediated contribution or found it 
to be negligible, we show that taking universal supersymmetry breaking masses 
at the Planck scale, rather than at the gauge unification scale as is customary, 
leads to the gluino contribution being more significant and in fact sometimes 
even larger than the chargino mediated contributions when fj, > and tan (3 
is of order 1. The impact is greatest felt when the gluinos are relatively light. 
Taking the universal boundary condition at the Planck scale also has an effect 
on the chargino contribution by increasing the effect of the wino and higgsino- 
wino mediated decays. The neutralino mediated contribution is found to be 
enhanced, but nevertheless it remains relatively insignificant. 



The flavor changing decay b — > 57 is often an important test of new physics 
because it is rapid enough to be experimentally observable although it appears first 
at the one loop level in the standard model (SM), thus allowing new physics to add 
sizeable corrections to it. For example, the decay is useful to limit parameter space in 
the minimal supersymmetric standard model (MSSM) |], H, H, d, §, H, 0, H H, [HI HI 



This is an especially useful tool if certain constraints have already been placed on the 
MSSM. Since the decay vanishes in the limit of unbroken supersymmetry, the relevant 
constraints pertain to the terms in the Lagrangian that softly break supersymmetry 
(SUSY). The general soft SUSY breaking scalar interactions for squarks and sleptons 
in the MSSM are of the following form: 

V soft = Q(Au\u)U c H 2 + Q(Ad\d)D c H\ + L(A E X E )E C H 1 + h.c. 

+ Q j m 2 Q Q + U c ^m 2 r U c + D cj m 2 D D c + L j m 2 L L + E c] m 2 E E c , (1) 

where Q and L are squark and slepton doublets and U c , D c , and E c are squark 
and slepton singlets. It is most frequently assumed that soft breaking operators are 
induced by supergravity, and that these operators have a universal form which is 
generation symmetric and CP conserving. Usually this assumption is coupled with 
that of grand unification to provide further constraints on the model's parameters. 
For purposes of simplicity the universal boundary condition is traditionally taken at 
the grand unification scale, even though the soft-breaking operators would be present 
up to near the Planck scale. The universal soft SUSY-breaking boundary condition 
is described by the following parameters: the scalar mass mo, the gaugino mass mi/2, 
and the trilinear and bilinear scalar coupling parameters A and B, respectively. 



Despite convention, some recent papers ]12|, [13], [14], [15[ have demonstrated im- 
portant implications of taking the boundary conditions at the Planck scale Mp and 
evolving them down to the scale Mq where grand unification is broken. The impor- 
tant difference between taking the universal boundary condition at the Planck scale 
versus at the GUT scale is due to that the top quark mass is known to be large 
(~ 174GeV) and grand unification causes some fields to feel the effects of the top 
coupling by unifying them into the same multiplet with the top. For example in 
SU(5) grand unification, Q, U c , and E c together transform as the 10-representation. 
Taking the universal boundary condition at the GUT scale is ignoring the fact that 
the soft-breaking parameters run above the GUT scale. One may at first think that 
any effects of grand unification would be supressed by powers of 1/Mq, but it has 
been demonstrated that such effects rather depend on In (M P /M G ) |16| . One surpris- 
ing result, which is given in refs. [13|, [15[, is that the predicted rate for the lepton 



flavor violating decay fi — > e7 may be only one order of magnitude beneath current 
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experimental limits. This is found even in SU(5) grand unification, where only neu- 
tralinos are available to mediate the decay at the one-loop order. Although due to the 
soft-breaking masses being assumed flavor blind and the unitarity of the CKM mixing 
matrix one might naively expect the partial width to vanish, this does not happen 
because operators with strength of the top coupling cause the third generation scalar 
fields contained in the 10- representation, in SU(5), to be considerably lighter than 
the corresponding fields of the other two generations [[Tj|. Below the grand unification 
scale, the only squarks effected by the top Yukawa coupling are the top singlet and 
the third generation SU(2) L squark doublet. 

We will apply these facts to demonstrate that in some of the same regions of 
parameter space where chargino corrections are important, one may also expect size- 
able corrections from gluinos. Previous calculations of b — ► sj in the MSSM have 
routinely either ignored the contributions mediated by gluinos || £|, ||, |6], [7j, [UJ or 
found them to be less important 0, H 0] than we will find them to be when we 
take the universal SUSY soft breaking boundary condition at the Planck scale rather 
than at the GUT breaking scale, as is cutomary. This is due to the fact, as explained 
above, that taking the universal boundary condition at the Planck scale leads to a 
different low energy SUSY spectrum than taking it at the GUT breaking scale does. 
For some parameter space, not only is the contribution from the gluinos important, 
but it is also comparable to that of the charginos. The fact that the top squark soft 
breaking masses are lighter when the universal boundary condition is taken at the 
Planck scale rather than the grand unification scale is, of course, likewise felt by the 
wino and higgsino-wino mediated decays. 



The standard model amplitude for b — > has been derived in refs. [ 171 , [TBI, and 
expressions for the additional MSSM amplitudes have been derived in refs. |], ||, [|, 
with ref. [p]] containing the first and most complete derivation. The leading-order 
QCD corrected partial width is given by its ratio to the inclusive semi-leptonic decay 
width in the following form |Tj| : 

E (6 -> 37) 6a \V t * s V tb \ 2 , 2 

7 = 7~T ... ,2 \ c 7 \ m b)\ , 2 

T{b-*cev) 7rg(z) \V cb \ 

where a is the electromagnetic coupling, g(z) = 1 — 8z 2 + 8z 6 — z 8 — 24z 4 lnz with 
z = m c /mb = 0.316 ± 0.013 is the phase-space factor, and the inclusive semi-leptonic 
branching ratio is BR(b — > cev) = 0.107. The ratio of the CKM matrix entries, for 
which we will use the experimental mid-value, is |V t *Vtf,| 2 / \V c b\ 2 = 0.95 ± 0.04. The 
QCD corrected amplitude cj(rrib) is given as 



c 7 (m 6 ) = r? 16 / 23 
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c 7 (M w ) - -c 8 (M w ) 



1 - rr 2/23 
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+ $> i 7 ? \ (3) 
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with dj and 6j being given in ref. r\ = a s (M w ) / a s (rrib) , for which we will use 

r] = 0.548. The terms ^{Mw) and c$(Mw) are respectively A 7 , the amplitude for b — > 
S7 evaluated at the scale M w and divided by the factor A® = 2Gp^Ja/8iT 3 V t *Vtbmb 
and A g , the amplitude for b — > divided by the factor A^a s Ja. The effective 
interactions for b —>■ sj and b —>■ sg are given by 

A 

L eff = (A-sa^PnbF^ + A^P^G^) + h.c. . (4) 

In calculating C7(M^y) and Cg(Mw), we will however use the conventional approxima- 
tion of taking the complete MSSM to be the correct effective field theory all the way 
from the scale Mq down to Mw, and hence ignore threshold corrections. 

As previously stated, normally A 1 is taken to be approximately the sum of A^ , 
A^ , and A* . In such a case, the charged Higgs contribution adds constructively to 
the SM amplitude. On the other hand, the chargino amplitude may combine either 
constructively or destructively with the other two, and in some cases may even cancel 
the charged Higgs amplitude. Even though the squarks strongly couple to the gluino, 
the contribution from the gluino mediated diagrams are considered negligible because 
the three generations of down squarks d iL belonging to Qu are conventionally assumed 
to have degenerate soft-breaking masses at the GUT breaking scale. However, the 
mass parameter mj~ is reduced by a small amount relative to mj~ for the first two 
generations in running the mass parameters down from the GUT scale, and b-squark 
mass matrix has off-diagonal entries proportional to m;, as given in the following 
equation: 

2 = ( m | L3 + m b ~ \ ( 2M w + M z) cos 2/5 rn b (A b + fi tan (3) \ 

m ' b ~ \ 3 m b (A b + fitaa/3) mf + m\ + | (M w - M|) cos 2(3 ) ' U 

where tan/3 = t^/^i is the ratio of Higgs vacuum expectation values and /i is the 
coefficient of the Higgs superpotential interaction ilHxH?,- These two effects make 
the b-squark eigenvalues somewhat different from the down squark masses of the 
other two generations, however the total effect is insignificant compared to the mass 
splitting that takes place in the stop sector due to the size of the top quark mass. 
(See, for example, Fig. 8 in ref. pOfl .) For this reason, and because the chargino 
contribution includes an often highly significant higgsino mediated decay, the chargino 
contribution to b —>■ 57 is found to be very important for some regions of parameter 
space, while the gluino and neutralino contributions are conventionally either found 
or assumed to be of little signifigance when the universal boundary condition is taken 
at the scale Mq. 

Now, we will perform the calculation with the universal boundary condition taken 
at the Planck scale and run the soft breaking parameters from there down to the 
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weak scale. In the following discussion we will use the one loop renormalization 
group equations for the gauge couplings, top Yukawa coupling, and soft breaking 



masses (See refs. |12], ^1). We will also use the exact analytic solutions, in the 
form derived in ref. Jl5|], to these one loop equations. We will use the conventions 
for the sparticle mass matrices and the trilinear coupling parameter A4 as found in 



ref. [SUSY]. (A { -> -A, t in the RGEs and RGE solutions of ref. |5|.) We will use 
a s (Mz) = 0.12. For the purpose of illustration, we will consider the specific case 
where the Planck scale trilinear scalar coupling A = 0, tan/3 = 1.5, X t (M G ) = 1.4, 
and the grand unification model is the minimal SUSY SU(5) model. In another 
paper [p^] , we will look at larger regions of parameter space, including large tan/?, 
for SO (10) grand unification. If X t (Ma) is reduced significantly then so also would 
be the effects we are discussing. For our chosen values of a s (M z ), the M G scale top 
coupling, and tan /3, the top quark pole mass is 168 GeV and the SM branching ratio 
is about 3.0 • 10~ 4 . For larger tan/3 the gluino and neutralino contribution would 
be greatly increased 0, but at the same time this would tend to occur for sparticle 
masses where the chargino contribution to b — > 37 is large enough to rule out the 
region of parameter space. 

Since the off diagonal terms in the b-squark mass matrix are much smaller than 
the diagonal ones and give relatively only a small contribution to the mass splitting 
between the b-squark mass eigenstates, we choose for simplicity to take the b-squark 
mass eigenstates to be the soft breaking masses. (See ref. fll4"| .) The two types of 
diagrams that contribute to b — > 57 and b — > sg in SU(5) with dn running in the 
internal loop are shown in Fig. 1. The internal fermion line represents either a gaugino 
or a neutralino propagator. To derive the contributions to A 1 or A g , one must sum the 
graphs with an external photon or gluon, respectively, attached in all possible ways. 
It is possible and simplest to work in a basis in which Au, the soft breaking squark 
masses and the trilinear couplings A are always diagonal in generation space [[14] . 
The masses of the first two generations of squarks dn are essentially equal to their 
soft breaking masses, which receive renormalization effects only from gaugino loops, 
and are hence degenerate. The soft breaking mass of b^ is much smaller than that of 
the other two generations, as we will see, since the bi belongs to the same multiplet 
as the top above Mq. Noting that there is no mixing at the &#-6#-gaugino vertex 
introduced by SU(5) grand unification and using the unitarity of the CKM matrix 
V, one may express the contributions to A^ by the gluinos [l], |2|, |15| as follows: 

A\ = -C(R)e D M^{- V Q - J 3J y V Qli/ 3J 



7 v ' u w a 2 l Ml 



3 

2 



+ Vb 1 ( A d + /itan/3) [G (mf^, m| 3 J - G (mf^ 



bf?' Qil 
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+ ^ 1 (A b -A d )G(ml R ,ml 3L )}, (6) 
where G is given by 

G(m 2 ,mi) = — ^4 ■> , (7) 

V 1 2 l M z m\-m\ w 

with 

gi(r) = — — -[2 + 3r — 6r 2 + r 3 + 6r lnr] , (8) 

6(r — l) 4 

<fe(r) = ~* [r 2 -l-2rlnr], (9) 

where we have used = —1/3, C(R) = 4/3, and rj b = m b (m b ) / m b (M w ) , which 
we take to be rj b = 1.5. The analogous expression for neutralinos may be obtained 
from the above expressions by working in the neutralino mass eigenbasis and noting 
that the first term in Eq. (6) comes from the diagram of Fig. la with the bino 
propagator and the other terms come from Fig. lb with the bino-wino propagator. 
(See, for example, ref. [15fl.) We have not included the diagram with the higgsino- 



wino propagator. Because neutralinos do not interact with gluons, one finds that 
simply A*° = A*° /eo- On the other hand because a gluon can attach to the gluino 
propagator, one finds 

4 = - h!,g 2 ^h 2 } + ^-(2- C{G)/C{R))A^ , (10) 

with C(G)=3 and 

h^r) = (l-6r + 3r 2 + 2r 3 -6r 2 lnr)/6(r-l) 4 , (11) 
h 2 (r) = -(-1 +4r - 3r 2 + 2r 2 lnr)/2(r - l) 3 . (12) 

Notice that the gluino contribution to b — > sg can be highly significant. 

The mass of Qsl may be expressed in terms of the first generation squark mass 
m% as 

QlL 

,2 2 ( T Iz 



where I G = (3/8vr 2 ) jjg A 2 (m 2 H + 2m 2 0a + A 2 ) din//, and I z /2 is the analogous con- 
tribution obtained from running the scale down from Mq to My/. The integrals Iq 
and Iz may be obtained from the analytic one loop solutions in terms of m 2 ^ and 
Mz(Mw) as follows: 



0.71M 3 2 , I z « 0.19m| iL + 0.17M 3 2 , (14) 
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where we have taken Mp = 2.4 • 10 18 GeV. One can also find that the relevant weak 
scale trilinear scale couplings are Ab ~ — 1.38M 3 and Ad « — 1.55M 3 . We calculate 
the parameter /j, at the tree level and find /i 2 « 1.0m| - 0.038M| - 4200GeV 2 . 

To illustrate the relative sizes of the seperate contributions to the b — > 57 ampli- 
tude, we plot 

r Ai ee A;/A™ (15) 

versus gluino mass for different values of ^q 1l in Fig. 2 for the case fi > 0. Notice 
that the gluino contribution to Ay is sometimes comparable to that of the charginos. 
This happens when the gluino mass is light. Note that the neutralino contribution is 
insignificant as usual. In Fig. 3, we plot the resulting branching ratio for b — ► 57. The 
dotted line corresponds to the calculation neglecting gluino and neutralino mediated 
contributions, while the solid line represents the full calulation. In both cases, we have 
included the SM, charged Higgs and chargino corrections as found in || , which work 
very well for low tan j3. When the glunino mass is 150 GeV, the gluino contribution 
can increase the branching ratio by as much as about 20-percent. This happens 
when there is strong destructive interference between the charged Higgs and chargino 
mediated amplitudes, and the gluino mediated amplitude then contributes to the 
branching ratio mainly through its cross term with the SM amplitude. In Fig. 4 and 
5, we show the analogous situation with the universal boundary condition taken at 
the GUT scale. Notice that when gluino masses are light, the gluino contributions 
are about one-third the size as when the boundary condition is taken at the Planck 
scale. The effect of taking the universal boundary condition at the Planck scale has 
only a small effect on the total chargino contribution for the parameter space shown 
here, however we find that for other nearby regions, for example with tan/3 = 2, the 
effect of making the chargino contribution more positive, but smaller in magnitude, 
is more appearant. In Fig. 6a and 6b, we plot the branching ratio as a function 
of the Mq scale gaugino mass M 5 g for curves of constant m for the cases where 
the universal boundary condition is taken at the Planck scale and at the scale Mq, 
respectively. Notice that for the curves with uiq > 0, the branching ratios for the 
complete calculation in the two cases differ by about 10-percent when M 5 g =60 GeV, 
which corresponds to a not very light gluino mass of about 170 GeV. When /i < 
and tan/3 is of order 1, we find the contribution to be much less important due to a 
strong destructive interference between the two diagrams in Fig. 1. 

In conclusion, if one is to calculate the decay rate for the flavor changing process 
b — > 57 in a SUSY GUT with SUSY breaking communicated by gravity above the 
GUT breaking scale in the form of soft breaking mass terms, it is essential to include 
the GUT scale renormalization group effects. An important result of including these 
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renormalization effects is that the gluino contribution to the decay rate can now no 
longer be neglected when the glunino mass is relatively light. 

We thank N.G. Deshpande, E. Ma, S. Nandi, and A. Soni. This was work was 
supported by Department of Energy grants DE-FG03-94ER 40837 and DE-FG02- 
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Figure captions 



Fig. 1 : The two types of diagrams with dn running in the internal loop that can con- 
tribute signifigantly to b R — > s^7 and b R — > s^g. One must sum the graphs 
with an external photon or gluon attached in all possible ways. 

Fig. 2 : Plots of = A^/A^ versus gluino mass for different values of rriQ 1L for the 
case \i > and tan/5 = 1.5 with the universal boundary condition taken at 
the scale M P . The curves correspond to squark masses raj =200, 300, 400, 
and 500 GeV. The gluino masses for each curve range from 150 GeV to the 
corresponding value of mj^. For example, raj =200 GeV corresponds to the 
curve for which the gluino mass ranges from 120 Gev to 200 GeV. Figs. 2a, 2b 
, 2c, and 2d correspond to the charged Higgs, chargino, gluino, and neutralino 
contributions, respectively. 

Fig. 3 : Plots of the branching ratio of b — > 57 for the case of Fig. 2. The solid lines 
represent the calculation including SM, charged Higgs, chargino, gluino, and 
neutralino contributions. The dashed lines represent the calculation using only 
the SM, charged Higgs, and chargino contributions. The curves represent the 
same squark masses as have been used in Fig. 2. 

Fig. 4 : Same as Fig. 2, but with the universal boudary condition taken at the Mq 
scale. 

Fig. 5 : Same as Fig. 3, but with the universal boudary condition taken at the M G scale 
as in Fig. 4. 

Fig. 6 : Plots of the branching ratio for the case of \x > and tan/5 = 1.5 as a function 
of the Mq scale gaugino mass M 5 g for curves of constant m . Fig. 6a corre- 
sponds to the universal boundary condition taken at the Planck scale. Fig. 6b 
corresponds to the universal boundary condition taken at the GUT breaking 
scale. The solid lines represent the calculation including SM, charged Higgs, 
chargino, gluino, and neutralino contributions. The dashed lines represent the 
calculation using only the SM, charged Higgs, and chargino contributions. The 
curves represent, in descending order in the two plots, m =0, 250 GeV, and 
500 GeV. 
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